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X ABSTRACT 

Efficient  procedures  for  time-stepping  Galerkin  methods^for  approximating  the 
solution  of  a coupled  nonlinear  system  -for  c(X,T)~  and  p-  with  non- 

linear Neumann  boundary  conditions, of— ^he  form 

V • [a (x ,c)  {Vp-Y(x,c)Vd}]  X-V  • u - , x e fl  , t c (0,T] , 

V • Ib(x,c , Vp) Vc]  - u(x,c,Vp)Nv  Vc  = 'frwfr  + f2<c) » xefi  ' t e (0,T^  ' 


= g^Xjt) 

= g(x,t,c) 


c{x,0)  = cQ(x) 


, XtJB  , t f (0  ,T]  , 
, x e 3fi  , t £ (0 , T] , 
, x c £1  , 


where  £2  c 3R^,  are  presented  and  analyzed.  Systems  like  (*)  ar^possible  model 

systems  for  describing  the  miscible  displacement  of  one  incompressible  fluid  by 

A 

another  in  a porous  medium  when  flow  conditions  are  prescribed  on  the  boundary. 

The  procedures  involve  the  use  of  a preconditioned  iterative  method  for  approxi- 
mately solving  the  different  linear  systems  of  equations  arising  at  each  time  step 
in  a discrete-time  Galerkin  method.^  Improvements  in  starting  procedures  over  many 

methods  are  obtained.  Some  negativi^index  norm  results  are  obtained  which  allow 
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weaker  smoothness  assumptions  on  — than  in  some  previous  treatments.  Optimal 
order  convergence  rates  are  obtained  in  most  cases  for  the  methods  which  are  com- 
putationally more  efficient  than  standard  methods.  Work  estimates  of  almost  opti- 
mal order  are  obtained.  The  techniques  developed  are  also  applied  to  single  non- 
linear parabolic  equations  with  nonlinear  boundary  conditions. 

AMS(MOS)  Subject  Classifications:  65M15,  65N15,  65N30,  76.35 

Key  Words:  Galerkin  methods.  Error  estimates,  Nonlinear  boundary  conditions. 
Fluid  flow. 
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Significance  and  Explanation 

Many  physical  problems  in  which  fluid  flow  or  heat  flow  is  prescribed  across 
the  boundary  of  the  region  of  interest  can  be  modeled  by  nonlinear  parabolic 
partial  differential  equations  (or  systems  of  equations)  with  nonlinear  boundary 
conditions.  An  analysis  is  presented  for  a particular  model  system  which  has  been 
used  to  describe  the  miscible  displacement  of  oil  by  certain  chemicals  in  the 
chemical  flooding  of  oil  wells  to  produce  greater  recovery  from  reservoirs.  The 
techniques  extend  to  a wide  class  of  such  physical  problems. 

Some  results  for  single  parabolic  equations  are  extended  to  certain  coupled 
systems  of  partial  differential  equations  where  optimal  order  error  estimates  are 
obtained.  Next,  approximation  theory  results  are  obtained  which  allow  some  re- 
duced smoothness  assumptions  on  the  unknown  functions.  Also,  analysis  is  presented 
which  allows  starting  procedures  which  are  computationally  more  efficient  than 
in  many  earlier  results.  Finally,  significant  amounts  of  computation  are  saved 
by  an  iterative  time-stepping  procedure  for  approximating  the  solution  of  the  large 
systems  of  linear  equations  produced  by  a Galerkin-type  numerical  procedure. 

Instead  of  factoring  a different  large  matrix  at  each  time  step  to  solve  the  linear 
equations  exactly,  only  one  matrix  must  be  factored  and  used  as  a preconditioner 
in  an  iterative  procedure.  Very  few  iterations  are  then  required  at  each  time  step 
since  the  iterative  procedure  is  just  a stabilizing  process  for  the  underlying 
time-stepping  procedure.  A complete  error  analysis  is  presented.  The  paper  con- 
tains work  estimates  which  show  the  large  computational  savings  of  the  precondi- 
tioned iterative  technique.  Almost  optimal  order  work  estimates  are  obtained. 


Ihe  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC , and  not  with  the  author  of  this  report. 
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Sufficient  condition*  for  f2.14)  to  held 


pr»««kM  r*i  « ♦ ior.*  11  y pr«ror<4i»  lor.***  *>»*  -%*  ~pt  inq 


In  this  4*v#lof>  • priori  lr«M»  for  f*«  +rr  r%  ?-<  *rwj  V(p-p) 


thrn  there  nut  constant*  K - r <i, a,, a»,b,, ,1-0 6),  ror  ‘h«  fourth  tern  on  the  rtqht  of  (4.5),  we  obtain 


often  be  estimated  by  using  the  test  function  * - - Atd 
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2 y- 


(n-1) AtanV(p-p)n,Vcn.cn) | 


We  integrate  fay  parts  to  bound  T . As  In  (4.21),  *e  have 
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where  !1  c E , are  presented  and  analyzed.  Systems  like  (*)  are  possible  model 
systems  for  describing  the  miscible  displacement  of  one  incompressible  fluid  by 
another  in  a porous  medium  when  flow  conditions  are  prescribed  on  the  boundary. 
The  procedures  involve  the  use  of  a preconditioned  iterative  method  for  approxi- 
mately solving  the  different  linear  systems  of  equations  arising  at  each  time 
step  in  a discrete-time  Galerkin  method.  Improvements  in  starting  procedures 
over  many  methods  are  obtained.  Some  negative-index  norm  results  are  obtained 
which  allow  weaker  smoothness  assumptions  on  than  in  some  previous  treat- 

ments. Optimal  order  convergence  rates  are  obtained  in  most  cases  for  the 
methods  which  are  computationally  more  efficient  than  standard  methods.  Work 
estimates  of  almost  optimal  order  are  obtained.  The  techniques  developed  are 
also  applied  to  single  nonlinear  parabolic  equations  with  nonlinear  boundary 
conditions . 
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